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A parallel current can destabilize a single flux line (FL), or an array of FLs. We consider the 
effects of pinning by point impurities on this instability. The presence of impurities destroys the 
long-range order of a flux lattice, leading to the so called Bragg glass (BrG) phase. We first show that 
the long-range topological order of the BrG is also destroyed by a parallel current. Nonetheless, some 
degree of short-range order should remain , whose destruction by thermal and impurity fluctuations, 
as well as the current, is studied here. To this end, we employ a cage model for a single FL in the 
presence o f impurities and current, and study it analytically (by replica variational methods), 
and numerically (using a transfer matrix technique). The results are in good agreement, and in 
conjunction with a Lindemann criterion, provide the boundary in the magnetic field-temperature 
plane for destruction of short-range order. In all cases, we find that the addition of impurities or 
current (singly or in combination) leads to further increase in equilibrium FL fluctuations. Thus 
pinning to point impurities does not stabilize FLs in a parallel current j z , although the onset of this 
instability is much delayed due to large potential barriers that diverge as j^^. 
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I. INTRODUCTION AND SUMMARY 

There is currently much effort devoted to understand- 
ing the nature of vortex phases in high temperature su- 
perconductors Q| . The flux lines (FLs) in a clean material 
prefer to sit on the sites of a triangular lattice forming 
an Abrikosov solid for fields H > H Cl . However, this 
mean field phase diagram is modified by the combined 
effects of thermal fluctuations and various types of disor- 
der, such as oxygen vacancies, columnar defects, or twin 
boundaries. Understanding the precise nature of these 
new phases is very important for the determination of the 
transport properties of these materials. In this paper, we 
focus on the combined effects of point-like disorder, such 
as oxygen vacancies, and a parallel current. (There are 
many studies that consider other types of disorder |^,|| ) . 

It is known that pinning to even a small amount of 
random impurities is sufficient to destroy the long-range 
translational order of the Abrikosov lattice |Q. If ran- 
domness is strong, the underlying order is destroyed at 
every scale, but for weak enough randomness, some order 
is expected to survive at short distances. Indeed, it has 
been suggested that there may exist a stable, dislocation- 
free Bragg Glass (BrG) phase at low magnetic fields and 
temperatures in the presence of weak point impurities 
H[|. It is important, both for technical applications of 
high T c materials and from a theoretical point of view, 
to understand the stability of this new phase with re- 
spect to an external current. For example, one may wish 
to control the melting transition of the BrG phase by 
some external parameter such as current. The case of 



a current perpendicular to the vortex lines was studied 
in Ref. 7j, where it was shown that the driven lattice 
may still maintain a topologically ordered moving BrG 
phase. Here, we study the case in which current is ap- 
plied parallel to the FLs. The instability of the proposed 
equations for the BrG suggests that the long-range topo- 
logical order of this phase is most likely destroyed by a 
parallel current. Using a Lindemann criterion and a cage 
model, we then discuss the destruction of any remnant 
short range order. This is in the same spirit as a recent 
calculation pf]| , in which we used a Langevin equation 
to examine the dynamics and fluctuations of a FL in the 
presence of an alternating parallel current. The basic idea 
is that even in the absence of true long-range order due 
to current or impurities, there remains a great deal of 
translational order in the vortex system at short scales. 
A "phase diagram" is then constructed to indicate the 
apparent limits of this "ordered" region. 

The response of FLs to a parallel current j z was stud- 
ied extensively even prior to the high-T c era H. While, a 
straight FL is stable in the presence of j z , any fluctuation 
of the line causes helical instabilities at long wavelengths 
H . A lattice of FLs is similarly unstable in the presence 
of an infinitesimal j z , if it penetrates the whole sample. 
As shown by Brandt |9) , for a finite penetration depth of 
current, a finite j z ~ e 1 ' 4 , is needed to cause instability, 
where e is related to the shear modules Cqq. Brandt has 
suggested || that FLs can also be stabilized for infinitesi- 
mal j z (which decays exponentially into the sample) due 
to pinning by disorder. However, this study predates 
the introduction of the BrG phase, which is more or- 
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dered than a flux liquid. We may thus ask if the BrG 
phase is stable with respect to j z l More generally, what 
is the behavior of the flux line in response to the cur- 
rent j z when point impurities are present? Since, the 
melting of the BrG phase is thought to be first order, 
the possibility of diverging longitudinal and perpendic- 
ular penetration depths of the current at the transition 
point is excluded. Nonetheless, we shall assume a cur- 
rent j z which uniformly penetrates the whole sample, 
and calculate the mean-squared fluctuations of the FL 
based on this assumption pl] , ^0| . This quantity is then 
used to examine the extent of short-range order, and to 
construct an apparent "phase diagram" in the magnetic 
field-temperature (B — T) plane. 

Two models used for the study of FLs with disorder 
are introduced in Sec. [il]. The first model describes the 
elastic behavior of the vortex lattice in the presence of 
impurities, and is appropriate for the study of its long- 
wavelength fluctuations. The true long-range order of the 
lattice is destroyed by such fluctuations, and possibly re- 
placed by the weaker topological order of the proposed 
Bragg glass . The second model considers a single FL 
confined by a harmonic cage potential that represents in- 
teractions with the other FLs. This model provides an 
approximate measure of the extent of short-range fluc- 
tuations, and in conjunction with a Lindemann criterion 
can be used to construct a phase diagram §. The addi- 
tion of cu rren t to these models is straightforward. 

In Sec. Ill replica variational methods are used to es- 



timate FL fluctuations in the presence of impurities and 
current. The BrG phase is described by a variational 
ansatz that breaks full replica symmetry. We propose a 
natural extension of this solution in the presence of cur- 
rent. This solution, however, has unbounded fluctuations 
and suffers from the same long- wavelength instabilities as 
the uniform lattice. We interpret this as a signal of the 
destruction of BrG order by the current. (Naturally, we 
cannot rule out other potential variational solutions with 
other forms of remnant long-range order.) However, it is 
quite likely that some degree of short-range translational 
order survives for weak currents and impurities. Using 
the cage model, we then estimate the fluctuations of a sin- 
gle FL. To this end, we generalize the replica variational 
ansatz proposed by Goldschmidt jl2| to include a cur- 
rent. The mean-squared fluctuations of the FL are found 
to grow with increasing current and/or impurities. Con- 
sequently, in the absence of the cage potential, a single 
FL is unstable in the presence of an infinitesimal current, 
and is not stabilized by the impurities alone. 

We test the predictions of the replica calculations (for 
the single FL in a cage) using a transfer matrix method 

For a given realization of 



which is introduced in Sec. IV 



randomness, this method provides numerically exact re- 
sults for the extent of fluctuations of a discretized FL in 
a cage potential. While it is difficult to provide an exact 
correspondence with the corresponding continuum limit, 
the results are qualitatively in agreement with those of 
the replica variational ansatz. In particular, we again 



find that both impurities and currents enhance transverse 
fluctuations, and that the FL is unstable in the absence 
of the cage potential. 

The conclusion that pinned FLs are unstable, even to 
infinitesimal current j z , may appear as somewhat coun- 
terintuitive: The force due to the current on a portion 
of the FL is then much less than the attractive poten- 
tial from impurities. The resolution lies in the size of 
the energy barriers that separate metastable pinned con- 
figurations. In Sec. [v| we argue that the size of these 
barriers diverges as j where the exponent /i is related 
to collective pinning effects. At finite temperatures, an 
exponentially long time is required for activated hopping 
over these barriers. Thus the onset of the instability can 
be much delayed due to pinning by impurities. 



II. MODEL HAMILTONIANS 

A. Flux lattice with impurities 

To discuss the large-scale deformations of the Bragg 
glass, we employ the elastic description introduced by 
Giamarchi and Le Doussal ||. If the displacements of 
the FLs with respect to their equilibrium positions, Ri, 
are denoted by the two-component displacement field 
u(Ri,z), then the energy cost of deforming the lattice, 
in the continuum limit, is given by 



It. 



2^ J (2tt) 3 



u a (q)$ a p(q)up(-q) 



(1) 



+ / tPxV(x)p(x), 



where a, (3 — 1,2 label the coordinates. The first term 
in the above equation denotes the elastic energy, and 
the g-integration is over the Brillouin zone. The sym- 
metries of the lattice constrain the form of the elastic 
kernel & a f3, for example resulting in two elastic mod- 
uli as q — ► in an isotropic system. For the sake of 
brevity, following Ref. we shall assume the simpli- 
fied form of $ Q/ g(<7) = cq 2 S a p. The second term in 
Eqn. (Q) is related to the potential produced by the point 
impurities. In the limit where many weak impurities 
act collectively on a FL, the resulting disorder is mod- 
eled by a Gaussian random potential with correlations 
V(x)V(x') — A(|f — r'\)S(z — z'). Here, r and z denote 
coordinates perpendicular and parallel to the FL respec- 
tively, while A(r) is a short range function of the scale 
of the superconducting coherence length. The density of 
the vortices at a given point x is also given by 



p(x) = ^ § (r- Ri - u(Ri,z) 



(2) 



The average over the quenched disorder can be per- 
formed using the replica method. After averaging over 
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disorder, the n-th moment of the partition function at 
a temperature T is the same as that of the replicated 
Hamiltonian pi, 



Ao 
2T 



+ E^ cos { /? -[" a ( a; )~^( a; )]} 



(3) 



where a,b — 1, • • • , n refer to the replica indices. Due to 
the lattice periodicity, averaging over disorder introduces 
both harmonic (Ao) and anharmonic (A#-) interactions 
between the replicated fields, where K is any reciprocal 
lattice vector, and Ak is the Fourier component of A(r). 

Generalizing the above Hamiltonians to include a cur- 
rent j z , parallel to the magnetic field, is straightforward, 
and achieved by adding the term <fioj z u x d z il/2. In gen- 
eral the current j z may be inhomogeneous, decaying from 
the surface into the bulk. To further simplify the calcu- 
lations, we shall assume that the current uniformly pen- 
etrates the sample |pj]| . This assumption enables us to 
write the replicated Hamiltonian as 



Heff — - 



<1 



(2tt) 3 



J2 G o^p(l)<(lH(-q) 



Ao Q a a b 

7^d a u a dpup 



(4) 



where G ^(g) = cq 2 S a s + i<t>oj z q z e a g. Pro perties of this 



Hamiltonian will be discussed in Sec. Ill A 



B. The Cage Model for a single FL 

To examine fluctuations of a representative FL at short 
scales, we shall use the cage model. The collective inter- 
actions of the other FLs is then approximated by a simple 
harmonic potential. The applicability of the cage model 
to the BrG is discussed in Refs. |3 10 1 . We again consider 
the limit where many weak impurities act collectively on 
the FL, so that the disorder can be modeled by a Gaus- 
sian random potential. Thus, the Hamiltonian for a sin- 
gle FL, whose position is denoted by the two-component 
vector f(z), is 



H 



dz 



£l 
2 



+ 7T r 2 + V 0, r) + 



t>a]z 
2 



x <9 Z 



(5) 



Here, E\ — £o/7 2 is the line tension of the FL, j 2 = 
m z /mx 3> 1 is the mass anisotropy, £0 = (0o/47rA) 2 , 
and 4>q = hc/2e. The spacing of the FLs is indicated by 
a, while A is the penetration depth 0. The spring con- 
stant is approximated in the dense limit of B > H Cl , 
by fi w £0 1 a 2 ■ The last term of the Hamiltonian is gen- 
erated by the Lorentz force due to the parallel current. 
The random potential is assumed to have the following 
correlations [p||l2]], 



(V(z, r)V(z', f ')) = Ae^ 3 4 2) (f - f')S(z - z'), (6) 
where 



f')~l/(2^)exp -(f-f')72£ 



(7) 



The dimcnsionlcss parameter A is a measure of the 
strength of the disorder, and £ is the vortex core diameter 
(coherence length) . 

Goldschmidt p2[ has used replica and variational 
methods to compare and contrast the effects of point and 
columnar disorder in the above Hamiltonian. In Sec. Ill B 



we shall generalize his results to the case of a non-zero 
current. 



III. THE REPLICA VARIATIONAL METHOD 

A. Instability of the Bragg Glass 

In the absence of disorder, the elastic Hamiltonian 
of Eqn. (0) has normal mode energies of cq 2 ± 4>oj z q z . 
Clearly, there are negative energies at small q z for any 
finite j z , indicating the instability of the lattice to long 
wavelength deformations. We ask if a similar instability 
is present in the presence of disorder, when the vortex 
lattice is replaced by a Bragg glass. To this end, we 
shall first briefly review the approach of Giamarchi and 
Le Doussal ||, and then generalize it to a finite j z . 
It is convenient to first decompose the propagator 



0,aP 



G oXp = ( C( 7 2 + Mzqz)[(S a p + ie a p)/2] (8) 
+ {cq 2 - (j> j z q z )[{S a p - ie a p)/2] 



U r aP T u r aPi 



where P^g R are projection operators (to left and right 

circular polarizations), which satisfy PapPp^ = SijS a -y. 

The essence of the variational approach is to find an 
optimal Gaussian trial Hamiltonian of the form 



Ho — 



d 3 



G- L p, ab (q)u a a (q)u (-q), 



0) 



2(2tt) 3 

^ere G~l b = Gj' L (q)P^ + G^^P* The 



generic form chosen for the trial propagators is G ab 1,1 = 
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SabGg 1 ' 1 — o l ah . The matrices o-^[ R act as effective caging 
potentials due to disorder, and may in principle stabi- 
lize the system. The best variational approximant to the 
Hamiltonian of Eqn. (^) is obtained by minimizing the 
variational free energy F = F + (H e f / — Ho), and results 



L/R 



A* 



^ 2T K 6XP 

K 



K 2 



where 



B 



L/R 



T 



d 3 q 

(2tt) 3 



G L b l%) 



(10) 



(11) 



These equations have been discussed extensively in 
Rcf. p|. Remarkably, they have solutions which break 
the replica symmetry. The full breaking of replica sym- 
metry is best represented by replacing the symbols ab 
with a continuous variable v £ (0, 1). The previous equa- 
tions for matrix elements can now be written as coupled 
integral equations 



K 



2T 



K 2 



(B L (v)+B R (v)) 



B L ' R (v) = 2T 



d 3 q 
(2^ 



G L / R (q)-G L / R (q,v) 



(12) 



(The propagators G a b(q) and G aa (q) are replaced by 
G(q, v) and G(q) respectively.) It is apparent from 
Eqs. @ and Q that <r L (v) = cr R (v), and we shall 
henceforth drop the indices L/R. A typical solution for 
cr(v) is a constant for v > v c , and an arbitrary function 
for v < v c , where v c is itself a variational parameter. The 
al geb raic rules for inversion of hierarchical matrices give 



B L / R (v) = B L ' R (v c ) 



+ / duj 



(13) 



d 3 q 



2Tct'(uj) 



where <j'(v) — dcr/dv, [cr](v) = vo~(v) — J Q V duja(u), and 
d 3 q 2T 



B L/R (v c ) 



{2n)3G -X, L/ R {q) + [(T]ivc) 



(14) 



To simplify calculations, we keep only the first term, 
with K = Kq, in the sum in Eqn. (pi]). (This is the 
single cosine model of the Ref. Q).) Eqs. (p"l|)-(p"4|) are 



now solved by first differentiating the expression for a(v) 
with respect to v. Eliminating o~'(v) from both sides of 
the equation, and performing the integrals over q, yields 



<j{v) 



Sue 3 ' 2 
K*T 



(15) 



for [a] > (jz/j o) 2 , and a(v) = for [a] < (j z /jo) 2 , where 
jo = y/ 4c/cf>Q. Taking one more derivative, and using 
[c]'(v) = va'(v), gives 



[*)= - 



(16) 



where v = K 2 T/ (Aire 3 / 2 ). 

Note that the above solution for [a](v) is always pos- 
itive, and hence incompatible with the initial condition 
of [a](v — 0) = 0. A potential resolution 1 is to assume 
that it goes to zero dis continuously for v less than some 



v c , i.e. 



, for v < v' c , 

H"') ! (i) 2 +(fe) 2 . iorv' c <v<v c , 
E , for v > v c . 



(17) 



Self-consistency of this solution with the condition 
[a](v) = va(v) - dujo-(u) requires v' c = v (j z /j ), 
and thus gives [a](v — v'+) = 2(j z /j ) 2 . We can as- 
sume that [er] is continuous at v c , leading to the condi- 
tion v c = vay^T. — (jz/jo) 2 - Another relation between 
these two unknowns is obtained by substituting v = v c 
in Eqn. (p^), leading to 



'E- 



fz Apg 4 
fo 



167TC 3 / 2 



exp 



M(B L (v c ) + B R (v c )) 



(18) 



where B L / R (v c ) is defined in Eqn. (fl4|), with [tr](u c ) = E. 
Solutions to these equations (for j z = 0) are extensively 
discussed in the Ref. ||. 

We shall now inquire if these solutions to the varia- 
tional equations are physically sensible. In particular, 
the variance of fluctuations at a wavenumber q is con- 
trolled by the correlation functions, 



G L l R (q) = G L a /R (q) 



1 dv 



[a](v) 



(G- 1 - L/R (q) + [a](v)) 
(19) 



Using the expression for [a] (v) , and performing the inte- 
grals over v, leads to 



1 Suggested to us by M. Mezard in private communications. 
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G L / R (q) = 



1 



cq 2 ± fojzQz 



Cl 



eg 2 ± <?W*gz + (j z /jo) 2 



c 2 (cq 2 ± ^ j z <7z 



(eg 2 ± + (j z /jo) 2 ) 3/2 

l-v c S 

w c eg 2 ± ^oj^fe + S 



(20) 



where ci = (l/u£ - l/v c )(j z /j ) 2 , and c 2 = 
[tan _1 (v c /i;£) — tt/4]/vo. For j z = 0, the leading be- 
havior as <f — > comes from the third term in Eqn. (20), 
and results in G L I R ~ 1 / j <gr| 3 . This is the characteristic 
signature of the BrG, quite distinct from the 1/q 2 scal- 
ing of an elastic material in the absence pinning, and 
leads to logarithmic growth of fluctuations in real space 
||. The modified singularity can be traced back to the 
v — > limit of the integral in Eqn. (19). A finite current 
removes this singularity as indicated by Eqn. (]l7|), and 
in the limit of q — > the integrals evaluated in the square 
brackets of Eqn. ( pp| ) go to a constant. Thus, apart from 
a multiplicative constant, the form of the correlator is 
preserved, and the solution suffers from the same insta- 
bilities as in the absence of pinning. (In other words, our 
Gaussian variational solution has some eigenmodes with 
negative eigenvalues.) We therefore conclude that this 
natural extension of the BrG in the presence of current 
is not stable. While it is possible that there are other 
variational solutions that do not suffer such instability, 
the most likely implication is that the (quasi-)long-range 
order of the BrG is destroyed by the parallel current. 



B. A single FL 

The mean-square fluctuations of each FL (due to short- 
wavelength fluctuations) are better estimated using the 
cage Hamiltonian of Eqn. (||). The average over the 
quenched random potential is again performed by intro- 
ducing n copies of the FL. After averaging over disorder, 
the replicated displacements (in Fourier space) of the FLs 
are coupled through the Hamiltonian 



dq z 



Once more, it is convenient to decompose the propagator 
in Eqn. ( |22|) into a sum of left and right circular polariza- 
tions as in Eqn. (||), where Gq 1,l / r = (eiq 2 ±4>ojzqz+^)- 
In terms of the projection operators, the inverse of the 
propagator for the trial Hamiltonian is simply 



Gab, a/3 — G a J b(q z )Paf3 + ^abilz) Pa/3 ■ 



(24) 



The trial propagator is written as G ab 1 a Aq z ) = 

SabG'^^qz) - a a b,oif3, with (Tab,af3 determined by min- 
imizing the variational free energy 



F = (H n - K) K - Tlog (j ' d[r a ]t 



-h n /T 



(25) 



Using Eqs. (|2l|), (|23|) and (|25|), and after some algebraic 
calculations one finds the free energy 



F 
2L 



constant + / |f £ G ^p(q z )G aa , af) (q z ) 

(x,f3,a 

| f ^trlnGfe) + i^/ p ($> abia J , (26) 

J ' a,b \ a / 



where the function f p is given by 
A^ 2 , 1 



fp(y) 



Air £ 2 + y ' 



(27) 



and 



B a b,af3 — T J [G aa ,af3 + Gbb,af3 — 2G a b, a p] ■ (28) 

The stationary solution for the free energy is obtained 
for <Tab,af3 = S a f3cr ab , with 



0"ah = Tpfv 



B. 



for 



b(^a) 



(29) 



Here, f p is the derivative of f p with respect to its argu- 



G^(?«)rS(«,)^(-9*) ( 21 ) ment, and^ /K are given by 



me 2 

47T 



exp 



a , 6 



2£ 2 



B 



L/R 



T 



dq z 
2vr 



G h ' H + G. 



L/R 
bb 



2G 



L/R 



(30) 



where 



(£i9 z 



■ fi) S a /3 + iq z <l>ojz£al3- 



(22) 



We now proceed to discuss the solution of the station- 
arity equations with the assumption of replica symmetry 
(RS), i.e. with a a b = & for ah off-diagonal elements. In 
the limit of n — > one finds, 



We again search for the best variational Gaussian Hamil- 
tonian h n , parameterized as 



G 



L/R 



G 



L/R 



5ab + o(G L J R )' 



(31) 



h n = 



dq z 
2(2?r) 



G- b ) a0 (q z )r a Mrl(-q z ). (23) 



Using Eqs. ( J30| ) and ([u]) and integrating over q z we ob- 
tain 



5 



ryL/R 



2T 



v^ 



m 1 



for 



Jz < Jc 



(32) 



and no solution for j z > j c , where j c = y/Aeifj,/^. (Note 
that in performing the integrals over q z from — oo to +00, 
we implicitly assume a lattice spacing smaller than the 
length V Jt/ei.) From Eqn. (p9|), the value of a is ob- 
tained as 



2 ~f< 



2T 



(33) 



It is interesting to note that in the limit of /1 — > 0, cr is 
equal to zero which means that there is no influence of 
disorder in the RS solutions in the absence of the cage 
potential. 

The transverse excursions of a long FL are given by 



(f 



?2\ 



T 



dg z 
(2n 



[G, 



L 
aa 



G. 



aa \ 



(34) 



Using Eqs. ( p4|) and (j3l]) and integrating over q z gives 



V^UVjcY 



(J 



2fi(l-(j z /j c f) 



(35) 



It is convenient to introduce dimensionless parameters 



T 

T = , and 

sot, 



B = —f-, (36) 



in terms of which we have /i = EqB/£? and a 2 = £, 2 /B. 
Using Eqs. (|3^) and (|36|), and rewriting Eqn. ( p5| ) in 
terms of these dimensionless parameters, we obtain 



« 2 VI (Jz/Jc) 2 

7 A [1 



(37) 



-3/2 



4nVB 



B(l-(j z /j c ) 2 ) 



-1/2 



Since in the real system we have B -C {jT) 2 , Eqn. j37| ) 
can be further simplified to 



jTVB 



7A 1 

4^ ( 7 f) 3 



(38) 



As expected, this relation predicts that the mean square 
displacements of the FL grow with increasing disorder 
and/or current. 

The mean square fluctuations from Eqn. ([38]) are plot- 
ted in Figs. 1(b) and 1(d), for zero and finite disorder, 
with the parameters 7 = 125, £ = 30A, £o£ = 1905°^, at 
T = 50°K, and with B = 250G. (Figures 1 (a) and 1(c), 
at the lower temperature of T = 30° K require the use 
of 1-stcp RSB, which will be discussed shortly.) These 



are the same set of parameters used by Goldschmidt in 
Ref. fL2fl . Indeed, we have essentially followed his ap- 
proach and notation throughout this section, generalizing 
to a finite j z . There is, however, one point of difference 
in calculating the extent of transverse fluctuations of the 
FL: We use Eqn. (HJ) which sums over contributions of 
all modes. Goldschmidt IQ considers fluctuations of a 
FL segment of length I* = V £ i/a*> which is similar to 
removing the contributions for q z < l/l* from Eqn. (|34|). 
Thus, for the same set of parameters, we obtain larger 
fluctuations (and a lower melting temperature). 



0.35 



0.25 



0.15 
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d 



0.8 . / . 
3 1 3c 



FIG. 1. The mean-square fluctuations of a FL with zero 
disorder (A = 0) for (a) T = 30° AT, (b) T = 50° K; and finite 
disorder (A/2tt = 0.2) for (c) T = 30° K, (d) T = 50° 

The next step is to look for variational solutions which 
break replica symmetry. The technical details of this 
calculation are presented in the appendix. Following 
Goldschmidt [p.2 , we find a solution with 1-step replica 
symmetry breaking (RSB) for temperatures below a j z - 
independent freezing transition point of 



T 



{£(> 



e/ 7 )(7A/2^) 1/3 . 



(39) 



The mean squared fluctuations with the 1-step RSB, 
again in the limit of B <C ("fT) 2 , are obtained from 




1 



1 



VcV 1 - Uz/jc) 2 

7 A 



1 



(40) 



1 



Vl + E/M - (Jz/Jc) 2 47T( 7 f )3 [1 - 0VJc) 2 ] 3/2 

1 n-vA 1 



Vi + s/m - Uz/j c y 



where the parameters v c and E are obtained by minimiz- 
ing the trial free energy expression given in the appendix. 
This minimization is carried out numerically, and the re- 
sulting points for T = 30°K < T cp « 35.4°^ are plotted 
in Fig. 1 (c) (and can be compared to those with zero 
disorder at T = 30°K in Fig. 1 (a)). 



G 



The calculated FL fluctuations are now used to con- 
struct an approximate 'phase boundary' for the destruc- 
tion of short-range translational order. Accordin g to the 
Lindcmann criterion, 'melting' occurs when \J (r 2 ) w 
clcl, where cl is a constant of the order of unity. From 
the replica symmetric solution of Eqn. (38), we obtain a 



melting field 




(41) 



In the limit of zero disorder, the above equation repro- 
duces the result B m = c 4 L (l - / / (jf m ) 2 of Ref. go), 
and agrees with the expected behavior of B ~ 1/T 2 for 
zero current. Using cl — 0.16, and the same parameters 
as in Fig. 1, the phase diagram is plotted in Fig. 2. For 
the portion of the melting line for T < T cp , we can no 
longer use the above analytical form, but must instead 
rely on the numerical results obtained using Eqn. (40). 
Interestingly, the melting field shows a maximum around 
the freezing temperature, even for zero current. Indeed, 
recent experimental results [|l4j show a reduction of the 
melting field for low temperatures. 
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FIG. 2. 'Phase diagram' of FLs for (a) A/2tt = and 
U/jc = 0, (b) A/2tt = and j,/j e = 0.5, (c) A/2n = 0.1 
and j z /j c = 0, (d) A/2tt = 0.1 and j z /j c = 0.5. 



IV. THE TRANSFER MATRIX METHOD 

With or without RSB, the variational method predicts 
that in the absence of the cage potential (fj, = 0), the FL 
is unstable in the presence of infinitesimal j z , and cannot 
be pinned by the impurities alone. We may still question 



the validity of this conclusion, given the known short- 
comings of the variational method. For example, in the 
absence of current, this method predicts that fluctuations 
of a FL grow with its length L as (r 2 ) ~ L. In fact, it is 
known that pinning to impurities results in (r 2 ') ~ L 2< *, 
with £ > 1/2 jnjyujj . We thus check the results using a 
numerically exact transfer matrix algorithm applied to a 
discrete model of the FL. 

Consider a FL pinned at one end to the point (0, 0, 0) 
on a cubic lattice, and denote the net Boltzmann 
weight of all paths terminating at the point (x,y,z) by 
W(x, y, z). We consider only paths that are directed 
along the z-direction, and hence this weight can be cal- 
culated recursively as 



W(z + 1, x, y) = exp 



W(z,x,y) + T x 



V{z,x,y) ^{x 2 +y 2 ) 



T 



2T 



(42) 



oVi*/i'cW{z, x+l,y) + e- yj */ j '°W{z, x-l,y) + 



~~W(z,x,y- 1) + e~* j * j °W(z,x,y+ 1) 



where j' c = 2T/(f> Q , and T = exp(-£i/2T). Note that at 
each step the path is either straight, or moves by only one 
step in the transverse direction. Such transverse steps ac- 
quire a weight r <C 1 from the FL elasticity, as well as a 
direction dependent contribution from the current [ jl7j . 
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FIG. 3. The mean-square fluctuations of a FL with F = 1, 
confined to a cage of n/T = 0.2, as a function of the reduced 
current jz/j' c , for (a) zero disorder; and (b) binary random 
potentials chosen from (0.2, 2). 

Starting with the initial values of W(0, X, y) = 6 0x So y , 
successive applications of this recursion relation provide 
the weights W(x,y,z) at subsequent layers, which can 
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then be used to calculate the (exact) degree of fluctua- 
tions Or 2 ) = (x 2 + y 2 ), for a particular realization of the 
random potential V(x,y,z). The results are then aver- 
aged over different realizations of the disorder potential. 
When the FL is confined to a cage (fj, ^ 0), Of 2 ) quickly 
saturates to a constant value. Fig. || shows the depen- 
dence of the saturation point on the external current. For 
this figure, the random potentials V/T were taken from 
a binary set (0.2,2), and the results averaged over 500 
realization of randomness. Since the microscopic models 
are different, it is hard to make a quantitative compari- 
son with the replica calculations. The numerical results 
are certainly in qualitative agreement with the replica 
estimates depicted in Fig. [I]. It is numerically difficult 
to probe currents close to the critical value, as the large 
transverse excursions require simulating large lattices. 




FIG. 4. Mean-square fluctuations of a FL in a 200 x 50 x 50 
lattice, with (a) V = 0,j z /j c = (b) V / 0,j z /j c = (c) 
V = 0, j z /jc = 0.08, and (d) V # 0,j z /j c = 0.08. 



found to be unstable, as depicted in Fig. f|(d). Indeed, 
we find that even an infinitesimal current can destabilize 
the FL, and that binary impurities only further increase 
transverse fluctuations. 

Since the use of binary randomness in directed poly- 
mer simulations may lead to anomalous results related to 
percolation fll8[| , we also repeated the simulations using 
a Gaussian distribution for the disorder potential V/T. 
Such impurities are more effective in pinning the FL, and 
appear to stabilize it at small scales. However, by vary- 
ing the width of the lattice or z, we see that the FL goes 
to the edges of the box, indicating its instability. Fig. || 
depicts the dependence of transverse fluctuations on the 
strength of Gaussian disorder. The fluctuations at short 
scales are governed by disorder, while at longer scales the 
current dominates, pushing the FL outward. 




FIG. 5. Transverse fluctuations for j z /jc ~ 0.1, with 500 
realizations of the random potential V/T taken from a Gaus- 
sian distribution of variance (a) 0, (b) 2, and (c) 4. 



The most interesting point is the behavior of the FL 
in the absence of the cage potential. As can be seen from 
the various limiting curves for fi = 0, plotted in Fig. ||, 
fluctuations no longer saturate, but continue to grow with 
z. When both the current and disorder are zero, the FL 
behaves as a directed random walk. The straight line in 
Fig. ||(a) corresponds to Of 2 ) cx t, as expected in this 
case. The addition of point disorder changes the scaling 
of transverse fluctuations Jig] . As indicated in Fig. ^(b), 
the average of Or 2 ) grows faster than t. The results for 
the binary set ({0.2, 2}) of random potentials V/T, aver- 
aged over 500 realization of randomness, are consistent 
with (r 2 ) cx t 2 ^ , with £ 0.6. This is in agreement with 
previous results for directed polymers in random media 
p5[ . As indicated in Fig. [|(c) the addition of current 
to the pure system (V = 0) leads to instability, and the 
transverse fluctuations become unbounded. When both 
binary disorder and current act on the FL, it is again 



V. ONSET OF THE INSTABILITY 

The results of replica and transfer matrix computa- 
tions can be summarized as follows: A parallel current j z 
causes instabilities in a flux lattice, or in a single (uncon- 
fined) FL. Point impurities appear to only increase the 
extent of fluctuations, and do not provide stability. On 
the other hand, if we start with a vortex system pinned 
by strong point impurities at low temperatures, it is diffi- 
cult to imagine how an infinitesimal j z can lead to drastic 
changes in configuration. In this section we discuss the 
length scales, barriers, and time scales that control the 
onset of the instability. 

Let us first consider a single FL in the absence of 
point impurities or a confining potential. A normal 
mode of wavenumber q z has energy proportional to 
G l ' L (q z ) = £i<zf i^oizfe) indicating a band of unsta- 
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ble modes for q z < <j>oj z /e\. The instability occurs only 
on segments of length longer than I* ~ £i/(<f>oj z )- (The 
divergence of £* as j z — ^ is cut-off by a finite confin- 
ing potential.) The simplest equation of motion for the 
FL relates its local velocity to the force obtained from 
variations of the Hamiltonian by df/dt = —p,[5H/8r], 
where jl is a mobility . The characteristic time scale 
for the normal modes is now given by T~ 1 ' L ^ R (q z ) = 

£iG 1,L ^ R (q z ). The time required for the growth of un- 
stable modes diverges as r* ~ A £ i/(^ojz) f° r jz —> 0. 
The FL thus appears stable at short times, or at high 
frequencies ui. 

The corresponding analysis for a flux lattice in the ab- 
sence of impurities is very similar. The (inverse) time 
scales for growth or decay of normal modes are related 
to the corresponding propagators G 1,L ' R m Eqn. (||) by 
a mobility factor. The unstable modes occur over an ap- 
proximately spherical region centered on q z < <f>oj z /2c 1 
and the time scale for the onset of instability scales as 

In the presence of point impurities, we claim that 
the instability still occurs for length scales larger than 
I* ~ (<t>ojz)~ 1 - The analytical support for this claim 
comes from the replica analysis: The best variational 
propagators (as for example in Eqn. (]2C|)) are propor- 
tional to the pure one with the same range of instabili- 
ties. There is also numerical support from the transfer 
matrix calculations: The comparison of Figures 4(c) and 
4(d) indicates that the onset of instability is not substan- 
tially modified in the presence of the disorder potential. 
The location of the instability (as determined from such 
curves for (f 2 )) is plotted in Fig. 6, and indeed scales as 
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In response to the current, the unstable FLs must 
change configuration at a length scale £*. This is not 
easily achieved in the presence of collective pinning by 
point impurities, and the flux system requires the assis- 
tance of thermal fluctuations to overcome the pinning 
barriers |l],|l9| . For a single FL the strength of the collec- 
tive barriers grows with the length scale as Et, oc £ s with 
9 w 0.3 (^0|. As j z — > the barriers at the length scale 
for the onset of instability thus diverge as Eb oc j J M with 
[i = 9. The probability that the flux system has sufficient 
thermal energy to overcome such barriers is proportional 
to the Boltzmann weight exp(— Eb/ksT). The waiting 
time for such rare events grows as r* ~ r exp[(jT/j 2 ) M ], 
where tq is a microscopic time scale, and jx is a charac- 
teristic current scaling as T^ 1 ^ 1 . 

The above arguments essentially follow those of Ref. 
|fl9f which considers the response of a collectively pinned 
flux system to a current perpendicular to the magnetic 
field. They are equally valid for a collection of flux lines 
in the vortex glass or Bragg glass phase. The exponent 
9 depends on the glass phase considered, its exact value 
is not as well known as for a single FL. Thus while the 
equilibrium state of a flux system is destabilized by an 
infinitesimal j z , the time for the onset of such instabil- 
ity diverges. The instability will not be observed if the 
system is probed at times shorter than r*, or perturbed 
by an AC current j z (u>) with frequency lu > lo* ~ 1/r*. 
(The latter is in qualitative agreement with the phase 
diagram presented in Ref. (HI.) 
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APPENDIX A: REPLICA SYMMETRY BROKEN 
SOLUTIONS FOR A SINGLE FL 

In order to find the replica symmetry breaking (RSB) 
solutions, we parameterize a a b by a function a(v) defined 
for < v < 1. The propagators G a b(q z ) and G aa (qz) are 
also replaced by G(q z ,v) and G(q z ), respectively. The 
stationarity equations are now written as 



\{B L {v) + B R {v)) 



where 



Fig. 6: The behavior of the position for the onset of insta- 
bility as a function of the current. 



B L / R (v) = 2Tj \G L ' R {q z ) G L / R (q z ,v) 
Inversion of the hierarchical matrices gives , 



(Al) 



(A2) 
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G L l R {q z )-G L / R {q z ,v) 
1 



(A3) 



1 du 



v (G^' L/R (q z ) + Jv ^ G^ L/R (q z ) + [(t](w)' 

where [<r](w) = wcr(w) — J Q duj(j(uj) as before. Integration 
over g z now leads to 



7.9 / ^s(MM) 



(A4) 



where cj([(t]) = T/^ei (m+MWoW- 

For a solution with 1-step RSB, we assume P3lp3 



a(v) 



(To f° r w < v c 
G\ for u > v c . 



(A5) 



This yields [<j](v) = for w < v c , and = E = 

w c (ci — do) for i > > u c . We can substitute the above 
form i n E qn. (A3) and use the result to evaluate B L I R in 
Eqn. (A2). Depending on whether v is bigger or smaller 
than v c , we obtain two different values for B L I R , which 
when placed in Eqn. (Al), lead to the self-consistency 
conditions 



1* 

T" 



2T 1 2T(1 - u c ) 1 
2T 



si 



(A6) 



if 



where si — y/4ei(n + E) — <^oiI an< ^ s 2 — 
(for j z < y/4£\(ij, + E)/^>o)- One can evaluate the free 
energy as a function of the parameters E and w c of the 
1-step RSB solution, as 



(A7) 
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Tl-t, c l 






E 




— = const. - 
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2L 


4 v c 
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2 ii 



E 

1 + - 
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1/2 | 


/' 







2T 



1 - «c ; 
2T p 

1 — v r 



y/l + E/n-fJft 



Vcy/l-jl/jl 



y/l + Y./VL-jlljl 



where t = T / y/JHi. Variation with respect to E repro- 
duces th e co ndition E = v c (ai — (To), uq and o\ defined 
in Eqn. ( |A6p . Minimizing the free energy with respect to 
the v c , after some algebraic manipulations, leads to 



T 



T, 



(A8) 



Cp 



with T cp — (A/27T7 2 ) 1 / 3 . Since v c < 1, this solution ex- 
ists only for temperatures lower than a freezing temper- 
ature, which is independent of t he current and [i. With 
this value of v c inserted in Eqs. ( |A6| ), we can numerically 
find E for any temperature and current, and evaluate the 
FL fluctuations in Eqn. (40). 



The natural question is whether there are also solu- 
tions wi th 2 -step, or even full RSB. Taking a derivative 
of Eqn (A4) with respect to a, and using the functional 
form for / p , leads after further algebraic manipulations, 
to a linear relation between u{v) and [(t](i>). This is in- 
consistent with the relation [c](u) = vo~(v) — J dwa{u>), 
unless a is a constant (since er'(u) = is always a solu- 
tion). One the other hand, one can also show by taking 
two derivatives with respect to v, that the only possible 
solution is v = v c — T/T cp , consistent with the 1-step 
RSB result. We thus conclude that, at least for this par- 
ticular form for the disorder correlator (and hence / p ), 
there is no solution with continuous RSB. 
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